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ABSTRACT 

In this second paper we present the first Nbody cosmological simulations of strongly cou¬ 
pled Dark Energy models (SCDEW), a class of models that alleviates theoretical issues related 
to the nature of dark energy. SCDEW models assume a strong coupling between Dark Energy 
(DE) and an ancillary Cold Dark Matter (CDM) component together with the presence of an 
uncoupled Warm Dark Matter component. The strong coupling between CDM and DE allows 
us to preserve small scale fluctuations even if the warm particle is quite light (« 100 eV). 
Our large scale simulations show that, for 10 11 < M/Mq < 10 14 , SCDEW haloes exhibit 
a number density and distribution similar to a standard Lambda Cold Dark Matter (ACDM) 
model, even though they have lower concentration parameters. High resolution simulation of 
a galactic halo (M ~ 10 12 M e ) shows ~ 60% less substructures than its ACDM counterpart, 
but the same cuspy density profile. On the scale of galactic satellites (M ~ 1O 9 M0) SCDEW 
haloes dramatically differ from ACDM . Due to the high thermal velocities of the WDM 
component they are almost devoid of any substructures and present strongly cored dark mat¬ 
ter density profiles. These density cores extend for several hundreds of parsecs, in very good 
agreement with Milky Way satellites observations. Strongly coupled models, thanks to their 
ability to match observations on both large and small scales might represent a valid alternative 
to a simple ACDM model. 
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1 INTRODUCTION 

A large variety of data on both small and large scales points to 
the existence of two dark components in our Universe: dark mat¬ 
ter, which dominates the gravitational budget of collapsed objects 
like galaxies and clusters, and dark energy, which regulates the ex¬ 
pansion of the Universe at later times (Tegmark et al. 2006, Planck 
collaboration 2014). While there is a general agreement on the pres¬ 
ence of these two components, their nature is still very unclear and 
several possibilities have been put forward in the last years. 

The most common model for Dark Matter (DM) and Dark 
Energy (DE) is the Lambda Cold Dark Matter (ACDM ). In this 


* E-mail: maccio@mpia.de 


model DE behaves like a cosmological constant, being often in¬ 
terpreted as vacuum energy, and is therefore constant in space and 
time. Furthermore, DM particle velocities are assumed to be neg¬ 
ligible, as though they decoupled quite early being already non- 
relativistic. 

The ACDM model is quite successful in reproducing data on 
large and intermediate scales (e.g. Springel et al. 2005). On the 
other hand it still raises several questions like the fine tuning prob¬ 
lem, as the DE energy density is ~ (10~ 3O m p ) 4 (m p : Planck 
mass), and the coincidence or “Why now?” problem: why DE, neg¬ 
ligible during all cosmic history, became significant just now and 
late enough to allow non-linear structures to develop. 

In the attempt to alleviate these problems, the option of DE 
being a self-interacting scalar field with a tracker potential was 
suggested (e.g., Ellis et al, Wetterich 1988, Ratra & Peebles 1988, 
Brax & Martin 1999, 2000, see Ratra & Peebles 2003 for a review). 
Furthermore the possibility of a coupling between DM and DE was 
deepened by several authors (e.g. Wetterich 1995, Amendola 2000 
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Amendola & Tocchini-Valentini 2001, see Amendola et al. 2013 
and references therein for a comprehensive review). 

While tracker potentials are devoid of finely tuned scales, the 
coupling allows for a non-negligible fraction of DE at early times 
(e.g. Amendola 2000). As expected, DM-DE coupling modifies the 
equation of motion of a test particle in an expanding universe with 
respect to the pure Newtonian case (Maccio et al. 2004). 

In the past years several works studied the effect of a DM-DE 
coupling on structure formation (e.g. Maccid et al. 2004, Caldera- 
Cabral et al. 2009, Baldi et al. 2010, Li & Barrow 2011. Baldi 
2012). After Planck data release, Xia (2014) outlined that the ten¬ 
sion between CMB and Hubble telescope Ho estimates was elim¬ 
inated by the coupling option, finding a coupling constant /3 = 
0.078 ± 0.022. Other observables, however, are not significantly 
modified by such coupling and this makes it hard to disentangle 
coupled dark energy models from a more simple ACDM one (e.g. 
Baldi et al. 2012, Pace et al. 2015). 

The ACDM model also faces some tension with data on the 
scales of low mass galaxies where the central cuspy dark matter 
distribution predicted by CDM (e.g Navarro, Frenk & White 1997) 
seems to be a very poor match of the observed central cored dark 
matter distribution in dwarf galaxies (e.g. Moore 1994, de Blok et 
al. 2001, Kuzio de Naray et al. 2006, Oh et al. 2011, Salucci et al. 
2012 ). 

In the last years, there has been a mounting evidence that in¬ 
cluding the effect of a dissipational baryonic component, strongly 
modifies the central DM distribution (Governato et al. 2010, 
Maccio et al. 2012a) reconciling observations and simulations on 
dwarf galaxy scales (Di Cintio et al. 2014a). 

On the other hand these simulations have also shown that 
baryons are effective in altering the DM density profile only up 
to a certain mass scale (Governato et al. 2012, Di Cintio et al. 
2014b, Arraki et al. 2014, Onorbe et al. 2015) and that very low 
mass galaxies are expected to retain a cuspy profile. 

Unfortunately this is in contrast with recent observations of 
the DM distribution in the MW satellites, which seem to suggest the 
presence of cored profiles even at these very low mass scales, where 
the effect of baryons should be minimal. The presence of such 
cores might indicate that DM is warm (e.g. Dalcanton & Hogan 
2001, Colin et al. 2008). Several works have dealt with the effect 
of Warm Dark Matter (WDM) on halo structure (Bode et al. 2001, 
Avila-Reese et al. 2001, Knebe et al. 2002, Tikhonov et al. 2009, 
Schneider et al. 2012). Constraints on the mass of a possible WDM 
candidate can be obtained from the analysis of the matter power 
spectrum from the Lyman-a forest (Seljak et al. 2006, Viel et al. 
2005,2008, Boyarsky et al. 2009), current measurements suggest 
the WDM mass (for a pure thermal candidate) to be around (or 
above) 3.5 keV (Viel et al. 2013, Polisensky & Ricotti 2014). 

A WDM particle of 3^1 keV, however, yields no improvement 
with respect to CDM on small scales (Schneider et al. 2014). The 
situation is particularly hopeless for halo density profiles, which 
require a WDM particle of ~ 100 eV in order to create apprecia¬ 
ble central cores as the ones observed in the Milky Way satellites 
(Maccio et al. 2012b, Shao et al. 2013). 

In the light of these considerations it is worth to explore alter¬ 
natives to ACDM , aiming to a model which can deal, at the same 
time, with both theoretical and observational issues. 

In two recent papers (Bonometto et al. 2012, Bonometto & 
Mainini 2014) a new model was indeed proposed, which combines 
DM-DE coupling and Warm Dark Matter trying to overcome the 
problems of ACDM , both on the theoretical and observational 
side. Observational DM, in this model, is warm and light. An aux¬ 


iliary CDM component, coupled to DE, is however added, which 
never exceeds some permils of the total density, but naturally ex¬ 
erts a number of key effects though cosmic history. 

In the companion paper (Bonometto, Mainini, Maccio 2015, 
Paper I hereafter) we addressed the linear behavior of such mod¬ 
els also including a new option, that CDM-DE coupling fades at 
low z. In these models WDM is made of extremely low mass parti¬ 
cles, down to few tens of eV, which might, in turn, have significant 
effects on the central DM distribution in collapsed objects. 

In this second paper, we aim to present the first Nbody simula¬ 
tions of strongly coupled warm+cold models (SCDEW hereafter). 
Our goal is to study these models in the highly non linear regime 
probed by structure formation, both on large (^Mpc) and small 
(~100 pc) scales. 

The paper is organized as follows: we will first summarize the 
main features of our novel cosmological model in section { 2 }. we 
will then introduce our numerical codes to generate initial condi¬ 
tions and evolve the simulations 0 . we will then present results 
for large scales simulations and for zoomed simulations of Milky 
Way-like objects and dwarf galaxies 0. we will then conclude with 
a discussion of our results and possible future developments 0 . 


2 THEORETICAL MODEL 

In this Section we briefly review the main features of SCDEW cos¬ 
mologies, more extensively discussed in the previous associated 
paper. Besides of baryonic and radiative components ( 7 ’s & i/’s), 
these models assume an ordinary uncoupled light WDM compo¬ 
nent, while DE is coupled to an ancillary CDM component which, 
actually, never plays the role that DM has in ACDM or similar mod¬ 
els. 

The starting point of SCDEW cosmologies is the finding that, 
during radiative expansion, an attractor solution exists, for coupled 
CDM and a scalar field components, keeping them a constant frac¬ 
tion of the cosmic density in primeval radiative eras. In the absence 
of coupling, these components would dilute oc a~ 3 and oc a~ 6 re¬ 
spectively. The time evolution of the different density parameters 
as a function of the expansion factor a in our SCDEW model is 
presented in figure |T] 

As shown in the Figure, the coupling, by allowing energy to 
flow from CDM to the field, puts CDM and DE on an attractor 
solution where both shall dilute oc a -4 , as radiative components 
do. For a coupling constant f3, the fair constant primeval density 
parameters of the two components, along the attractor, are 1/2/3 2 
and 1/4/3 2 , respectively. Values /3 ~ 10 are favored, so that the 
early contribution of CDM and field keeps steadily around or be¬ 
low « 0.2 %. While causing no harm to BBN or CMB data, this 
contribution therefore keeps non-negligible, possibly dating since 
the end of inflation. 

When this stationary solution is broken by WDM derelativisa- 
tion, the Universe naturally evolves towards the observed features. 
In particular, the field turns into the observed DE component. 

Accordingly, at variance from ACDM models, in SCDEW 
models DE has always been a small but non-negligible component, 
which eases fine tuning while the coincidence problem is also atten¬ 
uated, for a wide set of parameter choices, including those needed 
to obtain a reasonable fit to large scale data. The option of the field 
$ playing both the role of inflaton and DE is also not excluded. 

Coupled CDM however plays another key role. In the non- 
relativistic regime, as already known (Amendola 2000), its cou¬ 
pling to DE causes a strengthening of its self-gravity. This triggers 
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Figure 1. Time evolution of the density parameters (f 1i) of radiation, Dark 
Energy (4>), Cold Dark Matter, Warm Dark Matter and baryons in the 
SCDEW model. 


a mechanism allowing its fluctuations to restart WDM fluctuations, 
also on scales where they had formerly been erased by free stream¬ 
ing, as soon as WDM derelativizes. This revitalization of WDM 
fluctuations on scales below the free streaming one allows us to 
have very “hot” WDM, of the order of 100 eV. 

At more recent times (z <~ 50), besides of keeping a den¬ 
sity parameter < 10 - 1 -10 -2 of baryons, CDM could dynamically 
decouple from baryon and WDM fluctuations. 

This decoupling might arise naturally, however in our current 
parameterization we have imposed an ad-hoc decoupling at small 
redshifts, which causes no substantial change in the expected fluc¬ 
tuation spectra. Of course a full consistent model should provide a 
mechanism for the fading of the coupling at low redshift, which, 
for example, can be a conseguence of the dynamics of the scalar 
field. We refer the reader to Paper I (especially appendix A) for a 
thourough discussion of this issue while here we rather prefer to 
use a more phenomelogical approach. 

For the sake of definiteness, therefore, the linear spectra used 
to start the simulations in this work are those obtainable by setting a 
decoupling parameter d = 4 (see Paper I), yielding a full CDM-DE 
decoupling at 2 ~ 50, where simulations are started. 


3 NUMERICAL SIMULATIONS 

Numerical simulations have been performed with the Nbody code 
PKDGRAV (Stadel et al. 2001), while the initial conditions have 
been created with the GRAFIC2 code (Bertchinger 2001), which we 
have recently modified to allow a larger spectrum of cosmological 
models, as described in Penzo et al. (2014). The transfer function 
for the power spectrum have been produced with a modified ver¬ 
sion of CMBFAST that allows for an explicit coupling between the 
cold dark matter and the dark energy components (Bonometto & 
Mainini 2014). 


Figure 2. Linear power spectrum at z=0 for the ACDM model (blue) and 
the SCDEW model (red). The black thin line shows, for comparison, the 
power spectrum of a Lambda Warm Dark Matter model (LWDM) with a 
particle mass of 90 eV. 


The redshift-zero power spectrum for our specific SCDEW 
model is shown in figure [2] Despite the presence of a very warm 
component with a mass of 90 eV, the power spectrum does not dif¬ 
fer much from the expectation of a ACDM model with the same 
cosmological parameters. For the sake of comparison, in the same 
figure we also show the power spectrum of a pure WDM model 
for the same warm particle mass. As expected, it exhibit a dramatic 
reduction of the power on small scales. 

As is evident from Figure |T] at the initial redshift of our sim¬ 
ulations (z ~ 50), the CDM contribution to the gravitational po¬ 
tential is quite small. When combined with our assumption of a 
vanishing coupling at late times (see Paper I), this allows us to 
use the standard equations of motions for the Warm Dark Matter 
components and to neglect the presence of a Cold component in 
the simulation. On the other hand the initial power spectrum and 
the evolution of the background do account for the presence of all 
components, i.e., essentially, the effects of a fair time dependence 
of Dark Energy. 

Finally, in agreement with previous works on WDM (Colin 
et al. 2008, Maccio et al. 2012b, Shao et al. 2013), WDM particles 
were added the thermal velocity 


vo(z) 
1 + z 



/keV\ 3 

\m w ) 


kms 1 


(1) 


in agreement with Bode et al. (2001); here z is the redshift, h is 
the Hubble parameter in unit of 100 (km/s)Mpc _1 , fl w , m w and 
<)w are the density parameter, the mass and the number of spin 
states of WDM respectively. Their distribution function reads then 
( e Pw/Tw _j_ 1) —1 (p w and T w being the momentum and temperature 
of WDM), until gravitational clustering begins (Bode et al. 2001). 
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Table 1. Cosmological parameters for the ACDM and SCDEW models. 
Here 17$, Q c , Dw and Qf, are the density parameters of DE,CDM, WDM 
and baryons. All models share the same values for the Hubble parameter 
(h = 0.685), the spectral index ( n s = 0.968) and the power spectrum 
normalization (crs = 0.833). 


Label 


D c 


Df, 

m w [eV] 

SCDEW 

0.704 

0.001 

0.250 

0.045 

90 

ACDM 

0.704 

0.251 

- 

0.045 

- 


3.1 Large Box simulations 

We run two kinds of simulations, the first in a set of fairly large cos¬ 
mological boxes with the same resolution across the whole box. We 
have three different volumes with L = 20, 40, 90fi^ 1 Mpc aside, 
each of them containing 300 3 dark matter particles. For each box 
we run two simulations, one with our SCDEW model, and a refer¬ 
ence one done in standard ACDM for the same set of cosmologi¬ 
cal parameters (see table[T|for a complete list). We have used same 
values for the Hubble constant (h = 0.685), the spectral index 
{n s = 0.968) and the power spectrum normalization (crs = 0.833) 
in all models. 

In all simulations, dark matter haloes are identified using a 
spherical overdensity (SO) algorithm as described in Dutton & 
Maccio (2014). For our ACDM cosmology the virial density at red- 
shift zero is A(0) ~ 95.0 based on the fitting function of Mainini 
et al. (2003); we have used the same value for the SCDEW model 
as well. 

For each SO halo in our sample we evaluated the concentration 
parameter following the procedure outlined in Maccio et al. (2008). 
Briefly, we first determine the halo density profile, by using the 
most bound particle as the location of the halo center. We then com¬ 
pute the density (pi) in 50 equally spaced (in log) spherical shells. 
The minimum radius is the maximum between 1% of the virial ra¬ 
dius, .Rvir, or 3 times the softening length and the maximum ra¬ 
dius is 1.2f?vir- Errors on the density are estimated from the Pois¬ 
son noise due to the finite number of particles in each mass shell. 
We fit these density profiles with the NFW expression (Navarro 
et al. 1997): 

PnfwM = -j—. —(2) 
(r/r s )(l + r/r s ) 2 

here r„ is the scale radius of the halo, is normalization parame¬ 
ter, and per is the critical density of the Universe. Their values, and 
associated uncertainties, are obtained via a \ 2 minimization proce¬ 
dure using the Levenberg & Marquart method. We define the r.m.s. 
of the fit as 


77 ~ ln Pm) 2 , (3) 

i 

p m being the fitted NFW density distribution. Finally, following 
Maccio et al. (2007), we only selected relaxed haloes by requiring 
Prms < 0.5 and x 0 a < 0.07. Here p rms is the r.m.s. of the NFW fit 
to the density profile and a: 0 fF is the offset between the most bound 
particle and the center of mass, in units of the virial radius -Rvir- 

3.2 High resolution zoomed simulations 

We then selected a few haloes to be re-run at much higher resolu¬ 
tion, so allowing a more detailed study of the effect of our SCDEW 


Table 2. High resolution simulations properties. 


Halo 

Model 

M vir 

[h- 1 M 0 ] 

Rvir 

Softening 
|7i -1 kpc ] 

MW1 

SCDEW 

8.51 x 10 11 

1,573,013 

0.4 

MW1 

ACDM 

9.11 x 10 11 

1,685,767 

0.4 

D1 

SCDEW 

8.92 x 10 9 

1,501,683 

0.1 

D2 

SCDEW 

1.31 x 10 10 

2,205,387 

0.1 

D1L 

ACDM 

7.51 x 10 9 

1,264,309 

0.1 


models on the inner structure of haloes (e.g. density profiles) and 
on their satellite population. 

One halo comes from the 90 /i _1 Mpc box and has a mass 
similar to our own Milky-Way (« 10 12 /i -1 M 0 ). We will refer to it 
as MW 1. This halo was run with a mass resolution increase of 4096, 
reaching a mass per particle of 5.41 x 10 5 1 Mq and a softening 

of 0.4 /i -1 kpc . For this halo we performed three different runs, a 
ACDM one, two SCDEW, either with thermal velocities or without 
them. In order to identify bound subhaloes we used the AHlf] halo 
finder (Knollmann & Knebe 2009). 

Two further haloes (D1 and D2) have a much lower mass, 
similar to dwarf galaxies (« 5 x 10 9 h~ 1 Mg), and have been 
selected from the 20 /i _1 Mpc box. They were also zoomed in 
by a factor of 4096, obtaining in this case a mass per particle of 
5.94 x 10 3 ft -1 Mq and a softening of 0.1 /i _1 kpc . 

Unfortunately, no ACDM counterpart of these haloes could 
be run. As a matter of fact, on such small scales there are a lot 
of modifications to the tidal gravitation field when moving from 
ACDM to SCDEW, which alter the evolution pattern of the same 
initial Lagrangian region. In order to have a ACDM analogous of 
our dwarf SCDEW galaxies, we then turned to the 20 /i^ 1 Mpc box 
in ACDM and selected there a halo with mass and environment 
similar to the halo D2. We run it at high resolution, and will call it 
D1L. All the parameters of high resolution halos are summarized 
in tabled 


4 RESULTS 

We will start by presenting the analysis of large box simulations. 
Figure [3 shows the halo mass function in the two different models 
(obtained by combining all three different boxes of sizes 20, 40, 90 
/i _1 Mpc ). As expected from the Power Spectrum behavior (see 
figure 0. the two mass functions agree at high masses (large 
scales), while exhibiting some discrepancy as the halo mass de¬ 
creases. At intermediate masses the ACDM model (blue line) pro¬ 
duces more haloes with respect to the SCDEW one, while at the 
lowest masses probed by our simulations (M ~ 10 lo /i -1 Mg) 
there is an excess of haloes in the SCDEW model. Differences are 
small anywhere, never exceeding a factor of two. 

The concentration-mass plot shown in figure [4] exhibits more 
significant differences between the two models. This is expected, 
since the concentration parameter is more sensitive than the halo 
mass function to changes in the cosmological background (Maccio 
et al. 2008). While SCDEW shows concentrations similar to 
ACDM at the highest and lowest masses we probe, SCDEW ex¬ 
hibits significantly lower concentrations at intermediate masses: 

1 The Amiga Halo finder (AHF) can be freely downloaded from 
http://www.popia.ft.uam.es/AMIGA 
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Figure 3. Halo Mass function at redshift zero: ACDM is shown in red, 
while SCDEW is shown in blue. The lower panel shows the ratio between 
the SCDEW and the ACDM mass functions. 

5 x IO 10 — 5 x 10 13 h~ x M 0 . These lower concentrations are 
due to lack of power on those scales, in the SCDEW model; haloes 
then form at a lower redshift, so determining lower concentrations 
(e.g. Wechsler et al. 2002). 

As already noted in the Introduction, observations of low sur¬ 
face brightness galaxies and dwarf galaxies do suggest lower con¬ 
centrations with respect to the predictions of a standard ACDM 
model (e.g. de Blok et al. 2001; van den Bosch & Swaters 2001; de 
Blok & Bosma 2002; Swaters et al. 2003; Dutton et al. 2005; Gen¬ 
tile et al. 2005; Simon et al. 2005) and hence they favor SCDEW. 

4.1 Structure of Milky Way-like haloes in SCDEW 

We will now study in more details the dark matter distribution and 
the satellites population in our two Milky Way-like high resolution 
dark matter haloes. 

Figure[5]shows a logarithmic density map of the three versions 
of the MW1 halo: ACDM (left), SCDEW without initial thermal 
velocities (center) and SCDEW with thermal velocities (left). Even 
at naked eye, we notice a clear difference among the abundances of 
subhaloes in the three different cases. 

These differences are quantified in figure[ 6 ]where we show the 
subhalo mass function for the above three simulations. The lower 
number of satellites in SCDEW with respect to ACDM is due to 
the combination of two effects: (i) There is less power on the MW 
scale; this delays the formation of the halo and lowers the envi¬ 
ronment density at the formation time; in turn, this reduces the sub¬ 
haloes fraction (e.g. Maulbetsch et al. 2007). (ii) Furthermore, ther¬ 
mal velocities are high (WDM particles are 90 eV !) and this also 
suppresses the formation of low mass structures. 

A comparison of the three lines in figure [ 6 ] clearly shows that 
on the scale of the Milky Way, the former effect is more important, 
while thermal velocities only affect the tail of the subhalo mass 
function at masses M < 10 8 /i _1 Mg. 
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Figure 4. The concentration mass relation at z = 0. The blue points rep¬ 
resent ACDM , while the red ones show SCDEW. The black dotted line 
shows a simple lineal' fit (in log-log space) to the ACDM results. 

The subhalo mass function in the SCDEW model is quite dif¬ 
ferent from a pure Warm Dark Matter model, with the same WDM 
particle mass m v = 90eV. In fact, in the latter case, the formation 
of all substructures is suppressed, since the free streaming mass 
for such a model is of the order of few 10 14 /i -1 Mg (Zentner & 
Bullock 2003). On the contrary, in the SCDEW model, thanks to 
the coupling between CDM and DE, there are still enough sub¬ 
structures to host the formation of dwarf galaxies and, hence, not 
to violate constraints on the luminosity function of MW satellites 
(e.g. Maccio & Fontanot 2010). 

Figure[7]shows the radial density profile for our three realiza¬ 
tions of the MW 1 halo. The ACDM profile (blue) is slightly higher 
in normalization w.r.t SCDEW, while sharing a similar slope; the 
virial concentration parameters are 10.4 for ACDM and 7.1 for 
SCDEW; the difference is mainly due to the slightly larger virial 
radius in the ACDM run. The difference between the SCDEW runs 
with and without velocities is almost negligible. This is not surpris¬ 
ing given the halo mass and our choice of WDM mass (90 eV), 
we expect an effect on scales of few hundreds pc (e.g. Maccio 
et al. 2012), well below the resolution of our simulation. As we 
will see later, the situation is quite different for dwarf galaxies. 

4.2 Structure of dwarf haloes in SCDEW 

In order to investigate the effects of SCDEW on small spatial 
and mass scales we performed zoomed simulations of isolated 
dwarf haloes with masses of Mm = 8.92 x 10 9 (i _1 M 0 and 
Md 2 = 1-31 x 10 10 /i _1 M 0 . As already mentioned, due to the 
large differences in the tidal field between SCDEW and ACDM 
for such small scales, neither D1 nor D2 could be run in ACDM. 
We rather selected a different halo with a similar mass Mdil = 
7.51 x 10 9 /i _1 M 0 from the ACDM simulation, and run it at the 
same resolution; we dubbed it D1L. 

Figure[ 8 ]shows the density map of the D1 halo with and with- 
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Figure 5. Logarithmic density map of the MW1 halo in the different models, each panel is 400 h 1 kpc across, encompassing the halo virial radius. Left 
ACDM , center SCDEW without thermal velocities, right SCDEW with thermal velocities. 



Figure 6. Subhalo integral mass function for the MW1 halo. ACDM is 
shown in blue, while the solid and dotted lines represent SCDEW with and 
without thermal velocities respectively. 



Log (r/[kpc h" 1 ]) 


Figure 7. Density profile for the MW1 halo. ACDM is shown in blue, 
while the solid and dotted lines represent SCDEW with and without thermal 
velocities respectively. 


out thermal velocities. At variance from what happened on the MW 
scale, in the SCDEW run with thermal velocities all subhaloes have 
disappeared, and only the central halo has survived. This sets an 
effective free streaming mass ~ lO 7 ^^ 1 M@ for our model, equiv¬ 
alent to a warm dark matter (thermal) candidate of a few keV. 

The most interesting effect of the presence of a warm compo¬ 
nent can be seen on the density profile of our dwarf galaxies. Fig¬ 
ure [9] shows the density profile of the D1 (red) and the D1L (blue) 
haloes. The D1L profile exhibits the usual NFW-like shape with a 
cuspy central slope; the SCDEW profile without thermal velocities 
shows a NFW like behavior as well, with a cusp in the center, even 
if at a lower normalization. 

The SCDEW halo with thermal velocities, instead, shows a 
very clear cored profile. The size of the core is of the order of 
several hundreds of pc, in agreement with theoretical expecta¬ 


tions based on the conservation of phase space density (Maccio 
et al. 2012, Shao et al. 2013). Halo D2 exhibits a similar be¬ 
havior, as shown in figure [TO] The central slope of the profile 
is a = —0.25, substantially shallower than the NFW prediction 
(a = —1), even though the core is less pronounced. In this plot the 
blue line represents the theoretical expectation for an NFW pro¬ 
file in ACDM , assuming an average concentration for the halo D2 
mass. 

The flattening of the DM profile and the creation of a core 
on such small mass scales is one of the successes of the SCDEW 
model. Attempts of measuring the density profile of DM in the 
milky Way satellites suggests the presence of central cores (Walker 
& Penarrubia 2011, Amorisco & Evans 2012, Amorisco et al. 2013) 
even though there is still some debate on whether cores are re¬ 
ally needed to reproduce the kinematics of these satellites (Strigari 
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Figure 8. Density map, in logarithmic scale for the halo Dl, each panel is 80 h 1 kpc across, encompassing the halo virial radius. Left SCDEW with out 
thermal velocities, right SCDEW with thermal velocities 



Log (r/[kpc h- 1 ]) 


Figure 9. Density profile for the dwarf galaxy halo D1L (blue) and the Dl 
with velocities (solid red) and without them (dotted red). 


et al. 2014). If cores are indeed present in such small objects, it 
is quite unlikely that their creation can be ascribed to baryonic ef¬ 
fects (Garrison Kilmer et al. 2013), since the capability of baryons 
to modify the DM density profile is linked to the stellar/DM ratio, 
which is extremely low in these objects (Di Cintio et al. 14a). 

Any cosmological model able to create cores at this mass 
scales will definitely ease the tension with the observations and 
hence should be preferred to a simple ACDM model. 



Log (r/[kpc h" 1 ]) 


Figure 10. Density profile for the dwarf galaxy halo D2. SCDEW is shown 
as solid red, while SCDEW without velocities is shown by the dotted red 
line. The blue line is the ACDM theoretical expectation, for an NFW profile 
at this halo mass with an average concentration from figured 


5 CONCLUSIONS 

The ACDM model (based on DE with state parameter w = — 1 
and cold dark matter) is very successful in explaining the evolu¬ 
tion of our Universe and its redshift zero structure (e.g. Tegmark 
et al. 2007). 

Despite its success, this model faces problem on the theoret- 
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ical side, due to its inability to explain the very low value of DE 
density which, in turn, allows dark matter and dark energy densi¬ 
ties to be quite similar at the present day. Furthermore, such low 
value allowed DE to become relevant late enough, so to allow 
primeval fluctuations to turn into non-linear structures up to quite 
large scales. 

The ACDM scenario is also challenged by the possible pres¬ 
ence of cores in the observed dark matter distribution of Milky Way 
satellites (e.g. Walker & Penarrubia) since on such scales is very 
hard to invoke baryonic effects to reconcile observations and colli¬ 
sionless simulation predictions (Di Cintio et al. 2014a). 

It is then worth to explore alternatives to the simple ACDM 
model, possibly trying to find a model that can deal, at the same 
time, with both theoretical and observational issues. 

In this paper we presented a detailed non linear analysis of 
strongly coupled cold+warm cosmologies (SCDEW), while the lin¬ 
ear theory and the basic features of these models have been pre¬ 
sented in the companion paper (Bonometto, Mainini & Maccio 
2015, Paper I). 

Thanks to the coupling between Dark Energy and Cold Dark 
Matter, these models are able to ease the fine tuning and, possibly, 
the coincidence problem depending on the model parameters, prob¬ 
lems that plage a simple ACDM scenario. A very important feature 
of SCDEW is that the (uncoupled) warm dark matter component 
can be chosen with very low particle mass, as low as ~ 100 eV, 
since the strongly coupled CDM is able to regenerate fluctuations 
in the WDM component after its derelativisation (see Paper I). As 
a consequence, the power spectrum of the SCDEW cosmology is 
very similar to the one of standard ACDM, and hence in very good 
agreement with observations. 

We have run and analyzed several Nbody simulations per¬ 
formed in the SCDEW scenario and compared them with analo¬ 
gous simulations in the ACDM model. On large scales, SCDEW 
is quite similar to ACDM in halo distribution and number density; 
this implies that SCDEW shares the same success of ACDM on 
these scales. On the other hand, the distribution of mass in col¬ 
lapsed objects is less concentrated in SCDEW than in ACDM eas¬ 
ing possible tensions with observations (e.g. Salucci et al. 2012). 

By means of high resolution simulations, we then studied in 
detail the properties of haloes on the scales of our own Galaxy (the 
Milky Way) and dwarf galaxies. On MW scales, the SCDEW halo 
presents less substructures (30-40%) than its ACDM counterpart, 
possibly helping any “baryonic” solution (e.g. Benson et al. 2002, 
Maccio et al. 2010) of the so called missing satellites problem (e.g. 
Klypin et al. 2001). The total density profile is still well represented 
by an NFW fit, even though with a lower concentration parameter. 

Drastic improvements on ACDM are however found on the 
scale of dwarf galaxies Thanks to the large initial thermal velocities 
of the warm particles (due to the very low WDM particle mass of 
90 eV), the density profile in SCDEW haloes has a clear central 
core. The size of the core is of the order of several hundreds of 
parsecs, in good agreement with observations of the MW satellites 
like Sculptor and Fornax (e.g. Amorisco et al. 2013). 

SCDEW models are still in an infant state and more simula¬ 
tions and theoretical work is needed to fully test them at the same 
level as ACDM . For example we have negletted the possible for¬ 
mation of (very small) collapsed objects in CDM before the fading 
of the coupling. While these objects are expected to not sensible 
perturb the main gravitation potential, they could still leave some 
signature in the observational data (e.g. secondary CMB spectrum) 
and ought to be included in future numerical studies. 

Nevertheless our theoretical (Paper I) and Numerical (this pa¬ 


per) works clearly show that strongly coupled cold-warm cosmolo¬ 
gies are able to preserve the success of ACDM on large scales and 
to strongly improve the agreement with data on very small scales, 
simultaneously easing the theoretical short-comes of DE in ACDM, 
and thus appearing as a valid alternative to a simple ACDM sce¬ 
nario. 
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